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In our tunneling investigation using Andreev superconductor - normal metal - superconductor contacts 
on LiFeAs single crystals we observed two reproducible independent subharmonic gap structures at dynamic 
conductance characteristics. From these results, we can derive the energy of the large superconducting gap 
A L = (2.5 3.4) meV and the small gap A s = (0.9 -r- 1) meV at T = 4.2 K for the T l <? cal w (10.5 -r- 14) K 
(the contact area critical temperature which deviation causes the variation of Al). The BCS-ratio is found to 
be 2A_L/fcsTc = (4.6-5-5.6), whereas 2As jksTc <C 3.52 results from induced superconductivity in the bands 
with the small gap. 



The new class of superconducting rare-earth oxyp- 
nictides [T] is still not completely understood and there- 
fore requires further investigation. The layered LiFeAs 
(Ill-system) [2j is one of the few stoichiometric Fe- 
based pnictides which shows neither magnetic nor struc- 
tural transition but becomes superconducting at 18 K 
[21 S]. Band-structure calculations [51 [HI U\ show the 
Fermi surfaces for Ill-system to be comprised of quasi- 
two-dimensional (2D) hole cylinders centered at the T- 
point and electron ones at the M-point of the first Bril- 
louin zone that can be considered as two effective bands 
(so called minimal two-band model) [HI [5]. The total 
density of states at the Fermi level is formed mainly 
by Fe 3d-states pH [TTJ [H S3 E] ■ As was shown in 
[15] . the superconducting transition temperatures Tq 
for different types of iron-based superconductors corre- 
late with the total density of states at the Fermi level. 
This fact and the strong Fe isotope effect, which was 
reported by [TB] supports the phonon-mediated cou- 
pling importance [T7] in these compounds [TSl US] . The 
electron-phonon coupling is enhanced by an extended 
van Bove singularity [19) which was shown for iron pnic- 
tides and, in particular, for LiFeAs [3D]. In this work, we 
'present an investigation of the superconducting prop- 
erties of LiFeAs single crystals by means of Andreev 
spectroscopy of superconductor - normal metal - super- 
conductor (SNS) contacts, and the corresponding super- 
conducting gaps. 

The LiFeAs single crystals were obtained by self-flux 
method. The synthesis and investigation of the compo- 
sition and properties are detailed in |21j . A mixture 
of small lumps of the Li metal and powders of Fe and 
As in a molar ratio of Li : Fe : As = 3 : 2 : 3 was 
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placed into an alumina crucible. All work on the reac- 
tive mixture preparation was carried out in a dry box 
under argon atmosphere. The crucible was inserted into 
a niobium container, which was welded in argon at 1.5 
atm. The sealed Nb container was enclosed in a quartz 
ampoule. The sample was heated up to 1363 K, kept 
at this temperature and slowly cooled down to 873 K. 
At this temperature ampoule was extracted from the 
furnace and cooled in open air. The LiFeAs single crys- 
tals in the form of thin plates with lateral dimensions of 
(12 ± 6) x (12 ± 6) x (0.1 ± 0.05) mm 3 were separated 
from flux mechanically. According to the XRD, EDX, 
ICP MS, and nuclear quadruple resonance spectroscopy 
(NQR) obtained crystals have a stoichiometric compo- 
sition with a homogeneous distribution of elements over 
the entire sample |21) . 

Due to the strong hygroscopicity of LiFeAs, the 
rectangular plates were kept in glass capillaries. The 
"break-junction" technique [2H [23] was exploited to 
make SNS- Andreev contacts. The sample mounting on 
a spring holder was made in argon atmosphere in order 
to prevent the material from decomposition in air. The 
samples (thin plates with dimensions about (2 x 1 x 0.2) 
mm 3 ) were attached to the holder with two current and 
two potential leads by liquid In-Ga alloy. In order to 
measure on cryogenically clean surfaces, the holder was 
bended by a micrometric screw at T = 4.2 K causing a 
microcrack in the sample. The relatively small separa- 
tion distance between two superconducting banks pre- 
vents any extrinsic impurity penetration into the micro- 
crack. The current- voltage curves I(V), their deriva- 
tives dI(V)/dV and R(T) -dependences were measured 
by installation controlled by AT-MIO-16X (National In- 
struments) digital board [241 |2"5] . The dynamic con- 
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ductance spectra (dI(V)/dV) were obtained following a 
standard current modulation technique. 

Resistive measurements showed the superconducting 
transition of our Lii-^FeAs sing le crystals at T^ ulk rang- 
ing from 11 K to 17 K (see Fig.l), which may be caused 
by a minor affecting of water vapors and leads to a slight 
variation of lithium amount. As was shown theoretically 
in [55], while lithium deficiency turns Lii-^FeAs system 
into antiferromagnetic state (non superconducting) at 
8 = 0.5, structure of the Fe-As layers is not changed and 
c lattice parameter is decreased by ~ 14% [55]. Obvi- 
ously, the samples with critical temperature reduced up 
to ~ 0.7T™ Qa: (e.g. filled circles at Fig.l) are described 
by a small structure deviation from the stoichiometric 
one. To the best of our knowledge, a variety of physical 
measurements on the single crystals obtained from the 
same batch were carried out in [3] 14] I20 ] [2T ] [27 ] 128]. 
see also Table 1 concerning results on Tq and gap val- 
ues. The samples with the maximal T^ ulk w 17 K (see 
open circles at Fig.l, as well as data published in [2T]) 
demonstrated a relatively narrow superconducting tran- 
sition (1.2-5-1.5) K, which confirms high quality of single 
crystals. 

To evaluate the superconducting gaps two related 
methods were used: Andreev spectroscopy [5U] of 
individual superconductor-constriction-superconductor 
(ScS) 30 Sharvin-type contacts [3T] and intrinsic An- 
dreev spectroscopy (intrinsic multiple Andreev reflec- 
tions effect that usually exists due to steps-and-terraces 
presence at clean cryogenic clefts). Unlike SN- Andreev 
contacts, symmetrical ScS-contacts simplify the inter- 
pretation of the dynamic conductance spectra by us- 
ing theoretical model by Kiimmel et al. [30 . The 
main current-voltage characteristic (CVC) features of 
our ScS-contacts involve an excess current at low volt- 
ages and a subharmonic gap structure (SGS), showing 
series of dips of a dynamic conductance dl/dV at cer- 
tain bias voltages [301 [32] : 

V n = — , n = l, 2... (1) 

en 

due to the multiple Andreev reflections effect. In a 
case of a multi-gap superconductor several indepen- 
dent SGSs corresponding to the each gap should be 
observable (see Figs. 2-5). Such two distinct SGSs 
were reported earlier at Andreev spectra of Mg(Al)B2 
break-junctions [Ml 131 133, LaO .gF .iFeAs [SJ, 
GdO(F)FeAs [31 EZl EE], and FeSe [30]. The main 
advantage of a symmetrical SNS-contact is that at any 
T < Tc the gap energy can be defined directly from the 
bias voltages of the peculiarities at the dl /dV-cmve 
making use of the formula (1). The "break-junction" 



technique allows one to use such the advantages. The 
experimental CVCs obtained in this work are typical 
for the clean classical SNS-contact with excess-current 
characteristic J3H1 S2] , therefore, the model by Kiimmel 
et al. [30] is valid to describe our data. Strictly 
speaking, the sharpest SGS usually presents only 
at dI(V)/dV -characteristic of the most qualitative 
Andreev contacts with a small diameter a which is less 
than mean free path of the quasiparticles I (ballistic 
limit) 31 . In such a case one can observe several gap 
peculiarities, that increases the accuracy of the gap 
energies definition. Please note that if a « I, only the 
n = 1 minima would make a valuable contribution to 
the dynamic conductance spectra, which was the case 
for some LiFeAs Andreev contacts studied. 




T, K 

Fig.l. Superconducting transitions for LiFeAs single 
crystals measured before a microcrack forming. The 
dR(T) /dT-curve maxima were assig ned as T^ ulk . The 
bulk critical temperatures range from 11 K to 17 K 



A significant structure anisotropy of LiFeAs allows 
us to observe stack contacts (40j [41] representing a 
sequence of SNS-junctions and intrinsic multiple An- 
dreev reflections effect, which was observed earlier by 
the "break-junction" technique in other layered super- 
conductors (Bi-2201 [H], Mg(Al)B 2 [33 [32 I3S], and 
some Fe-based superconductors [Ml [3H [3H] along the c- 
direction). The bias voltage at ci/(l^)/dy-characteristic 
of such an array scales with the number of contacts N 
(in a stack) in comparison with dynamic conductance 
of an individual contact. Existence of natural stacks 
in c-direction on cryogenic clefts makes usage of such 
structures typical for "break-junction" technique. The 
array contacts provide reducing of surface defects in- 
fluence (which otherwise significantly affect the proper- 
ties of superconductor [33]) by N times. Therefore, the 
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magnitude of bulk gaps can be measured with a higher 
accuracy using stack contacts. 

Fig. 2 represents the excess-current CVC and the 
corresponding derivatives for two SNS-Andreev break- 
junctions (3 contact arrays normalized to a single junc- 
tion; were formed by a mechanical readjustment) in 
LiFeAs single crystal measured at T = 4.2 K. The main 
features at the dynamic conductance curve are located 
at Vi — 2A L /e pa 5 mV. The origin of the fine struc- 
ture is not yet clear and requires further investigation. 
The observed subharmonic gap structure (SGS) dips 
riL = (1,2,3) are seemingly related to the large gap 
(marked by black labels). Peculiarities at V pa 3.5 mV 
may appear due to an anisotropy of the Al order pa- 
rameter. However, we need further studies to confirm 
this unambiguously. The SGS bias voltages described 
by formula (1) (see inset of the Fig. 2) can easily yield 
the value A L w 2.5 meV (at the T l ° cal w 10.5 K) with- 
out any additional fittings. The "local" Tc is referred 
to as the Tc measured locally at the contact point with 
diameter a pa I and might differ from the Tc obtained 
by non-local methods. It is worth to note that location 
of Andreev reflexes were not moved under a mechanical 
readjustment of the contact while the contact resistance 
changed (see Ij^(V) and h{V) at Fig. 2). Thus, one can 
conclude that the gap values are still unchanged. This 
means high homogeneity of the sample superconducting 
properties in the measuring range. 

The dynamic conductance spectra of contacts j}ci5 
and jjd7 made on LFA6 sample (see Fig. 3) are pre- 
sented by light and dark-grey curves. Both contacts 
|d5 and jjd7 consist of two SNS-junctions in an array, 
while contact tfd7~ (plotted by black dashes) is a stack 
of three SNS-junctions (all curves are normalized to a 
single junction). The smearing of the gap peculiarities 
(Fig. 3) presumably originates from the comparatively 
large contact diameter (a pa I). Another smearing factor 
might be a temperature influence on the Fermi surface, 
which scales with ksT: for T = 4.2 K it yields an en- 
ergy uncertainty about 0.4 meV, which is of the same 
order of magnitude as the value of the small gap (see 
Fig. 4). These are, most likely, the reasons for the fuzzy 
small gap structure in the case. The complex structure 
of riL = 1 minima, as well as additional peculiarities 
in the range of V\ -r V3 voltages, can be caused by an 
anisotropy of the Al order parameter. By coinciding 
the 2Ai peculiarities at the stack conductance curves 
and normalizing them to a single junction (see Fig. 3), 
one can easily yield the large gap energy A^ pa 2.5 meV 

(jlocal _ 1Q 5 K )_ Thig yalue corresponds tQ all the CQn _ 

tacts of Fig. 3 and does not depend on the contacts resis- 
tance changing. Using our gap value, we derive a BCS- 




Fig.2. Normalized to a single junction CVC and dy- 
namic conductance of Andreev arrays Jjd7~ (bold light 
curve) and (|d8 (dark thin line) formed by 3 SNS- 
junctions in a stack (T = 4.2 K). The subharmonic 
gap structure (SGS) minima (til labels) define the 
large superconducting gap Al ~ 2.5 meV which is 
well-reproduced after the mechanical readjustment of 
the contact (the T l (? cal w 10.5 K). Inset: the linear 
SGS-minima bias voltages l/ n i(l/n)-dependence for the 
aforementioned characteristics 




Fig. 3. Normalized CVC and dynamic conductance 
(T = 4.2 K) of SNS-Andreev stacks of 2 junctions (con- 
tacts jjd5, ||ef7 plotted by light and dark solid lines, 
respectively) and of 3 junctions (contact jjd7~, black 
dashed line). The contacts were formed by a mechanical 
readjustment. The well-reproduced minima (til labels) 
define the A L » 2.5 meV (the Tg cal » 10.5 K) 



ratio for LiFeAs lA^lksTc ~ 5.5. In an attempt to 
consider those conductance spectra as 4 contacts (JJrf5, 
jjc£7) and 6 contacts (jJtfT^) in an array, one could get 
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2Ai/fcsTc 2.8 for the large gap that is inconsistent 
with the BCS coupling limit of 3.52. 




-8 -6 -4 -2 2 4 6 8 



Fig.4. Normalized to a single junction CVC and dy- 
namic conductance (T = 4.2 K) of jjc!4 (light curve) 
and j)d5 (dark curve) Andreev arrays (2 SNS-junctions 
in stacks). Two sets of Andreev minima (wl,s labels) 
define the large gap At 3.4 meV (the T^ cai « 14 K) 
and the small gap Ag = 0.9 -f- 1.0 meV. Linear back- 
ground was suppressed 

The SNS- Andreev spectra of stack contacts with the 
maximal T l ° cal « 14 K are shown at Fig. 4. The ar- 
rays representing two-contact stacks were formed by se- 
quent junction readjustment. The CVC and the dy- 
namic conductance comprising two sets of Andreev pe- 
culiarities in agreement with the formula (1) yield the 
values Al rs 3.4 meV and As w 0.9 meV (contact 
)|d4), As w 1 meV (contact j}d5). For the contact jjd5, 
the ns = 1 peculiarity was taken as the center of the 
doublet. Such a structure may be caused by a s-wave 
symmetry distortion of As due to the complex form 
of the hole-like Fermi surface sheets. The reproducibil- 
ity of the small gap minima (ns = 1) at dI(V) / dV- 
characteristics of Andreev stacks confirms the existence 
of the small superconducting gap and suggests it to be 
a bulk property of LiFcAs. The results presented at 
Figs. 2, 3 allow the same conclusion for the large gap. 
The minima ns = 2 for the small gap become observ- 
able on subtracting a strait line at the dI(V)/dV raise 
region near zero bias. 

The temperature affecting on Andreev peculiarities 
is shown at Fig.5a. The T l ° cal w 12.5 K is obtained 
from the linearization of the CVC derivative. The dy- 
namic conductance dips series (black labels) seemingly 
reveal the superconducting gap A rs 2.5 meV. It is easy 
to observe that A-minima marked as til = 1 start to 
close at T « 6 K. The gap temperature dependence 



A(T) (being in fact the Vi/2 positions plotted versus 
T) is presented at the Fig. 5a inset. The aforementioned 
T ~ 6 K appears to be approximately a half of the 
T local The A(T) was fitted by the single-band BCS- 
like curve and the two-band dependence (plotted in a 
framework of Moskalenko [33] and Suhl [33 equations 
for a superconductor with A^/As ~ 3.5 and a weak 
interband coupling). Both the theoretical curves tend 
to the T l (? cal = 12.5 K, in the same time, the best fitting 
is achieved if the two-band curve approaching T = K 
lies above the single-band one. This gives the average 
gap value about 2.5 meV. The fitting curves behavior 
indicates no evidence for an induced superconductivity 
in the condensate under consideration. Hence, this gap 
parameter is related to the bands with the large, "driv- 
ing" gap Ai and relevant neither to small gap, nor to a 
surface gap. 

The comparable temperature dependences of the 
large gap Al(T) for the contacts jjd5, LFA3 sample 
(see the dark curve at Fig.4) and j)c, LFA5 sample are 
shown at Fig. 5b. The experimental points were fitted 
by single-gap BCS-type functions. The reproducibility 
of the large gap behavior is obvious. Despite the differ- 
ence in the A^ and T(? cal values for the contacts, the 
corresponding BCS-ratios 2A^ /fceTc are close to each 
other and approximately equal to 5.6. The closer exam- 
ination brings out the Al{T) deviations from the BCS- 
like curves: all the experimental dependences slightly 
bend, which originates from nonzero interband inter- 
action with As-bands. As for the A$(T) temperature 
dependence represented by small circles at Fig. 5b (con- 
tact jjc£5, LFA3 sample), its behavior is typical for a 
"driven" gap. The similar phenomenon was observed in 
Mg(Al)B2 superconductor earlier [331 ES]- The differ- 
ence between Al(T) and As(T) courses becomes an- 
other evidence for existence of two nearly independent 
superconducting condensates described by the two gap 
values in LiFeAs. 

There are no published tunneling data (energetic pa- 
rameters obtained with a help of SIS-, NS-, NIS-, SNS- 
spectroscopies) on LiFeAs so far. It was revealed in 
the present work for the first time. Therefore, we com- 
pare our results with the data by the number of scien- 
tific groups in the Table 1. Angle-resolved photoemis- 
sion spectroscopy (ARPES) on LiFeAs single crystals 
[3J [13 [28], specific heat [571 ES], London penetration 
depth [2H1I3Z1I3H]: and microwave surface impedance [3S] 
yield similar results for the number and values of super- 
conducting gaps (see Table 1). The distinction in the 
BCS-ratio values should be taken into account. On cal- 
culating 2A/k B T c the Tg nset in [3j [2JJ [28j ES S3 SI ED] 
and T(V- lk in i4SJ have been used, which are average 
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Fig.5. a) Dynamic conductance of LiFeAs Andreev sin- 
gle SNS junction measured at the temperatures ranging 
from 4.2 K to T^° cal « 12.5 K. SGS defines the large 
gap Al ~ 2.5 meV (ul labels). The curves were shifted 
along the vertical scale for the sake of clarity. Inset: 
The large gap temperature dependence plotted using 
experimental data of Fig.5 (rhombs). The single-band 
BCS-like fit (dark line) and possible two-band fit (light 
line) were presented, b) The comparison of the large 
gap temperature dependences: sample LFA5, contact (jc 
(squares; T l c ocal = 12.5 ± IK, A L (4.2K) « 3meV), and 
sample LFA3, contact (td5 (circles; T^ cal = 14 ± IK, 
A L (4.2K) » 3.4meV, A S (4.2K) » lmeV; see Fig. 4). 
The point size reflects a gap value uncertainty. The 
single-band BCS-like dependences are plotted by dash- 
dotted lines 



values that may be higher or lower than the T(? cal in 
the point of gap defining due to a samples inhomo- 
geneity if any. By contrast, operating with the local 
Tc, one can calculate the BCS-ratio value more accu- 



rately. For the reason given above our BCS-ratio value 
exceeds that of other groups obtained by other meth- 
ods. Nevertheless, an excellent correspondence is re- 
vealed with the 2Ai/fcsTc value from [SO] (see Table 
1). For the large gap, our experimental data lead to 
ratio 2Ai,/A;bTc = (4.6 -j- 5.6) which exceeds the stan- 
dard one 3.52 for single-gap BCS superconductors in the 
weak coupling limit. At the same time, the BCS-ratio 
for the small gap 2As/ksTc <C 3.52 indicates the su- 
perconductivity induced by interband coupling with the 
"driving" bands (characterized by the large gap A^) in 
the bands with As due to k-space (internal) proximity 
effect [51 between two condensates, which resembles the 
situation in Mg(Al)B 2 [531 [33 [3S] , an d some Fe-based 
superconductors [Ml [36 [ |37 [ |38 [ [39 ] . 

In conclusion, it was found that I{V) and dI(V)/dV- 
characteristic of the SNS-Andreev break-junctions in 
LiFeAs single crystals with bulk critical temperature 
= (11 -J- 17) K could not be described in a frame- 
work of the standard single-gap model. Being observed 
by Andreev spectroscopy on LiFeAs for the first time, 
two independent subharmonic gap structures indicat- 
ing the existence of multi-gap superconductivity yield 
the following energies of two distinct superconducting 
gaps: A L = (2.5 3.4) meV and A s = (0.9 1) meV 
(at T = 4.2 K and the T l ° cal « (10.5 14) K). The 
uncertainty in the determination of the gap can be es- 
timated to be about 10% for large gap A^ and 20% 
for As owing to enormous raise of dynamical conduc- 
tance at low bias voltages. The estimated BCS-ratios 
2Ai/fcsTc = (4.6-T-5.6) exceeding the single-gap weak- 
coupling limit 3.52 can be caused by a strong coupling 
in the " driving" electron-bands [3] characterized by the 
large gap. On the other hand, the ratio for the small 
gap 2As/ksTc <C 3.52 hints on the importance of in- 
terband coupling (due to k-space proximity effect) that 
is to be taken into account when considering supercon- 
ducting properties of LiFeAs. 
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Table 1. Comparison of LiFeAs experimental data obtained by different techniques. *The mean Tc denned from 
the dR(T)/dT derivative was used. **The data range reflects an anisotropy of As order parameter [3]. 



Present data (SNS break- junction) 


Sample 


Contact 


T l ° ca \ K 


A L , meV 


Ag, meV 


2A L /k B T c 


2A s /fc s T c 


LFA6 


db, d7, d7~ , d8 


10.5 ± 0.5 


RJ 2.5 




5.5 




LFA3b 


c 


12.5 ± 0.5 


ss 2.5 




4.6 




LFA3b 


(24 


14 ±0.5 


ft 3.4 


ft 0.9 


5.6 


1.5 


LFA3b 


d5 


14 ±0.5 


RJ 3.4 


ft! 1 


5.6 


1.6 


Experimental data obtained on the single crystals from the same batch (as samples used by us) 


Ref. 


Technique 


^" set , K 


A L , meV 


Ag, meV 


2A L /k B T c 


2A s /k B T c 


Borisenko, et al. [3] 


ARPES 


18 


3.2 


1.5 -r 2.5" 


« 4.1 


1.9 4-3.2 


Stockert, ct al. [27] 


specific heat, ARPES 


17 


2.6 


1.2 


w 3.5 


ft 1.7 


Inosov, et al. |28| 


London penetration depth, ARPES 


17 


3.1 




4.1 




Experimental data published by other scientific groups 


Wei, et al. [36] 


specific heat 


17 


2.7 ±0.8 


0.5 ±0.2 


3.5 


1.2 


Song, et al. |37] 


London penetration depth 


17.5 


2.9 


ft! 1.3 


ft 3.8 


ft 1.7 


Kim, et al. g8] 


London penetration depth 


ft 17* 


ft 2.8 


ft 1.6 


« 3.8 


« 2.2 


Imai, et al. |49| 


microwave surface impedance 


17 


2.9 




4.0 




16.3 


3.0 


1.1 


4.2 


ft 1.6 


15.6 


3.0 


1.65 


4.4 


« 2.5 


Sasmal, et al. 1501 


vortex penetration, reversible magnetization 


14 -M5.3 


3.3 


0.6 


5.4 


1 



